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Abstract. We study directed transport in periodically forced scattering systems 
in the regime of fast and strong driving where the dynamics is mixed to chaotic 
and adiabatic approximations do not apply. The model employed is a square 
potential well undergoing lateral oscillations, alternatively as two- or single-parameter 
driving. Mechanisms of directed transport are analyzed in terms of asymmetric 
irregular scattering processes. Quantizing the system in the framework of Floquet 
scattering theory, we calculate directed currents on basis of transmission and reflection 
probabilities obtained by numerical wavepacket scattering. We observe classical as 
well as quantum transport beyond linear response, manifest in particular in a non-zero 
current for single-parameter driving where according to adiabatic theory, it should 
vanish identically. 
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1. Introduction 

The concept of pumping in nanosystems pLj has emerged from the paradigm of peristaltic 
pumps, simple idealized devices moving mass and/or charge in a defined direction 
as two or more parameters are cycled through a closed path [2]. As long as the 
driving is considered slow, adiabatic approximations apply. This enormously fruitful 
idea allows for a very general analytical treatment involving quantum scattering theory 
and Berry phases [HI HI El E] and has inspired a wide range of experimental realizations, 
including semiconductor nanostructures, e.g., quantum dots driven by time- dependent 
gate voltages [3 E] and more recently graphene sheets P [10]. Stimulated in turn by 
these developments, theoretical approaches have been extended towards quantum pumps 
with fast and strong driving, leaving the adiabatic regime — largely equivalent to linear 
response [6] — far behind [HI [121 (131 El [E], as well as towards including realistic many- 
body effects like dissipation and charge blocking. 

By contrast, aspects of complex dynamics involved in directed transport in strongly 
driven pumps have not yet been explored in depth. From the point of view of ballistic 
classical motion, peristaltic pumping amounts to a trivial dynamics of the transported 
particles, moving phase-space volumes through the device much like a viscous liquid 
[H ISl Uni |T^. As sufficiently coherent and controllable sources of ever faster drivings 
approaching the THz regime [H] become available, pumping in the realm of nonadiabatic 
and strongly nonlinear dynamics comes within reach. Scattering theory, already at the 
heart of quantum pumping, provides us with an appropriate theoretical framework which 
has developed far beyond the adiabatic regime: Irregular scattering [19] is a mature field 
with numerous applications from celestial mechanics [20] down to chemical reactions 
[2T| l22l [23] . comprising classical as well as quantum phenomena and approaches. 

With this work, we attempt a first survey into pumps operating in the regime 
of chaotic scattering and to give an overview of their principal features. As concerns 
implications for directed transport of the nonlinear dynamics induced by fast and strong 
forcing, there exists a pertinent case to be followed: The concept of chaotic ratchets [24] , 
in particular classical Hamiltonian ones amenable to direct quantization [251 126], has 
stimulated important insights into conditions and mechanisms of directed transport and 
its quantum manifestations in this regime. Similarly, based on a comparable body of 
results available on scattering at strongly driven potentials [271 l2H], we will point out 
peculiarities of pumping due to a complex dynamics, such as currents violating linear 
response ^9] and the possibility to achieve transport with a single driven parameter 
[131 [TSl l30] . Conversely, the field of irregular scattering gets enriched by considerations 
around symmetry breaking [31] and finds new applications with significant technological 
perspectives such as the generation of polarized (pure spin) currents [32] . Focusing on 
dynamical aspects, we however abstain from considering many-body effects in this work. 

In the first, classical part of the paper, we devise a family of elementary one- 
dimensional models, consisting of square potential wells which allow for two- as well 
as single-parameter driving, and analyze directed transport in these systems in terms 
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of phase-space structures. The quantum part briefly reviews Floquet scattering theory 
as basic theoretical tool and provides numerical evidence, obtained from wavepacket 
scattering, for directed currents largely owing to the complex underlying classical 
motion. The concluding section provides an outlook to a number of questions left 
open by our work. 



2. Classical chaotic pumps 



2.1. Models 

We seek models that (i) show irregular scattering, (ii) exhibit directed currents yet 
(iii) remain simple enough to permit an efficient numerical or even an analytical 
treatment, and (iv) are amenable to experimental realizations. Square potential wells 
and combinations of them almost ideally meet these conditions [271 i2Hl EH] : They allow 
for different types of external driving to induce chaotic behaviour and break symmetries 
as is necessary to achieve directed transport, and they closely resemble the potentials 
of semiconductor superlattices in the transverse direction. 




Figure 1. Models for single- and two-parameter driving. Panel a: Laterally driven 
square potential well as in eq. ([1]); b: position of the potential well (shaded) as a 
function of time for two-parameter driving, eq. ([2]); c: same for single-parameter 
driving, cqs. p|4|) . 



We consider a square well with constant depth but lateral driving, i.e., with the 
positions of its walls depending on time in some arbitrary periodic manner (fig. la). 



V{x,t) 



<S){x-x-{t)) - Q{x-x+{t)) 



(1) 



where Q{x) denotes the Heaviside step function and the wall positions x±{t) are periodic 
functions with the same period T. A lateral driving arises naturally from an AC voltage 
through a Kramers-Henneberger gauge transformation [331 [31]. A phase shift between 
the walls, equivalent to an oscillating width of the well, can be motivated, e.g., by a finite 
propagation time across the well of time- dependent perturbations or by a superposed 
gate voltage. Mass is set m = 1 throughout. 
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We assume either two independent harmonic driving forces ("two-parameter 
driving") with phase offset (fig. lb), 

x±{t) = + Xocos{ut±'ilj/2), uj = 2tt/T, (2) 

where a is the width of the well (in all that follows, a = 2), or alternatively ("single- 
parameter driving") a synchronous motion of both walls (fig. Ic) with the same periodic 
but otherwise arbitrary time dependence f{t), 

x±(t) = ±| + Xo/(t), fit + T) = fit). (3) 
2.2. Symmetries 

In periodically driven scattering, the phase of an incoming trajectory relative to 
the driving constitutes an additional scattering parameter JlS], besides the incoming 
momentum pin = \imt^_oomx(t). It can be defined, e.g., as 6* = wtin mod 27r, with 
iin the arrival time of the scattered particle, measured by extrapolating its asymptotic 
incoming trajectory till the origin x = ^29j. As 6 is typically beyond experimental 
control, directed transport is considered relevant only averaged over 6. In order to 
avoid a systematic cancellation due to counterpropagating trajectory pairs related by 
P: r — i- — r, r = {p,x), 6 ^ —9 (equivalent to w — > —u or t — )■ — t for harmonic drivings 
like ([2])), we have to break time- reversal invariance (TRI) of the potential. For eq. (|2]), 
this is achieved if 7^ 0, vr. In eq. (Ill), it requires a function f{t) without any symmetry 
or antisymmetry, f{t± — t) 7^ =t/(^) fo^^ reference times t±, which excludes in particular 
harmonic forces. In the sequel, we choose (cf. fig. [1]:) 

/(t) = cos(a;t) +7cos(2a;t-0), 77^0, 7^Omod7r/2. (4) 

For a two-parameter driving, even if) could be difficult to control. To make sure that 
currents do not even vanish upon averaging over if), we prevent the systematic pairing 
of trajectories related by P and ip — )■ —ip [29], e.g., by choosing functions for x±(t) not 
connected by any symmetry. For an asymmetric single-parameter driving such as (j4]), 
this problem does not arise in the first place. 

A relevant issue in the context of symmetry is the energy distribution of incoming 
and outgoing trajectories. In order to analyze deterministic transport processes based 
on irregular scattering, it would suffice to consider individual scattering events at a given 
energy, without assuming specific distributions in the asymptotic regions or reservoirs to 
keep them constant. The strongly inelastic scattering envisaged here, however, together 
with symmetry breaking suggest to have a closer look at this aspect. 

If, for example, we set = 7r/2mod7r ineq. (jl]) (such that f{7i/2ijj—t) = —f{'n/2uj+ 
t)), trajectories come in pairs T"{t") = Cr'(t'), C: {p,x,t) — j- {p,—x,—t) (fig. [2^). As a 
consequence, incoming and outgoing momenta are interchanged, ^('^(0^^) = Pout(in) ' 
the energies -Ein(out)- We illustrate the spread of Eout vs. Ei^ in fig. [2)3, c in terms of the 
probability density {d{Ei^, Eont))9- In panel (b), the trajectory pairing is reflected in a 
symmetry of the distribution with respect to the diagonal (the ridge along the diagonal 
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Figure 2. Symmetry and inelastic scattering. Panel a: Systematic pairing of 
transmitted (above) and reflected (below) trajectories for single-parameter driving 
with antisymmetric time dependence /(7r/2 — t) — — /(7r/2 + t), eqs. (|ll3l4p with 
(f) = 7r/2, see text and fig. [1]:. b,c: Probability density distributions ((i(i?in, £'out))e 
(color code) of outgoing vs. incoming energies for (b) antisymmetric driving as in (a) 
and (c) asymmetric driving, (j) — 7r/4. Other parameter values are Vb = 0.75, xq = 0.1, 
7 = 0.2, uj = 1.0. 



corresponds to elastic processes i^in = i?out)- It is absent (c) if the antisymmetry of the 
driving is broken choosing, e.g., = 7r/4 in eq. (jlj). 

Therefore, if i^in and i?out obey, say, Fermi statistics with the same chemical 
potential /i and < /i < E'^^^, the process i?in(= £^out) ~^ -^out(= -^in) is less probable 
than — )■ E'^^^^ and a symmetry-induced balance in the reflection probabilities (fig. [2^, 
lower panel) from either side is broken for inelastic processes, so that they now contribute 
to the current. More generally, inhomogeneities in the energy distribution, even if they 
occur identically on both sides, can well enable a finite contribution to transport of 
processes which otherwise would be suppressed by some symmetry. 

We show quantum results for Fermion reservoirs below (figlH]) but mainly consider 
transport, in classical as well as quantum calculations, as a function of E'in without 
averaging over this quantity. 

2.3. Phase-space structures 

In the context of Hamiltonian ratchets, the optimal condition to obtain directed currents 
is not hard chaos but a mixed dynamics with chaotic and regular regions coexisting in 
an intricately structured phase space [251 EE]. Pumps are no exception to this rule. 
We therefore focus on systems which pertain only marginally to the regime of chaotic 
scattering (defined through (i) the existence of a chaotic repeller in phase space, (ii) 
self-similar deflection functions, and (iii) exponential dwell-time distributions [T9]). 

In fig. |3] we show a set of Poincare surfaces of section, for two-parameter, eq. ([2]) 
(panels a,b) as well as for single-parameter driving, eqs. ( I3]l4l) (c,d). In all plots, we 
observe chaotic regions interspersed with regular islands, characteristic of the KAM 
scenario [35]. Typically, these islands surround periodic orbits bound in the scattering 
region and not accessible from outside (see, e.g., the inset in panel a). The phase-space 
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Figure 3. Poincare surfaces of section (PSS, defined hy x — 0, p > 0). a,b: two- 
parameter driving, eqs. (|1I2I) . c,d: single-parameter driving, eqs. (|1I3I4[) : a: with 
TRI, harmonic driving as in eq. ^ with -0 = 0, inset: periodic orbit underlying the 
prominent regular island, tracing the position x{t) (short dashed) between X-{t) and 
x+{t) (fuh lines); b,c,d: TRI broken with V = 7r/2 (b), 7 = 0.1 (c,d), (f) ^ (c), tt/2 
(d). Other parameter values as in fig. [2] Large void areas correspond to trajectories 
not contained in the initial ensemble or intersecting the PSS only a few times. 



structures clearly reflect the presence or absence of spatiotemporal symmetries of the 
underlying dynamics, as is evident comparing panels (a) (TR invariant driving) with 
(b),(c),(d) (TRI broken). 



2.4- Directed transport 

Moving stepwise from phase-space structures to currents, we analyze the relation 
between scattering and directed transport on two levels: (i) locally, the outcome of 
individual scattering processes in terms of asymmetries in reflection and transmission 
from either side, and (ii) globally, mean currents as functions of scattering parameters 
like i^in- They are obtained by averaging discrete outcomes (transmission vs. reflection) 
over, for example, the phase 9 (see subs. 12. 2p . 

{i)e = ^ do {[T-^{e) + R^^e)] - [T^~{e) + R-{e)]), (5) 

where T^^, R^"^, denote transmission and reflection probabilities, resp., between the 
asymptotic regions x — > ±00. Their familiar left-right symmetry is generally broken in 
the absence of a corresponding P-invariance of the Hamiltonian, so that, e.g., reflection 
from left to left can coexist with transmission from right to left at otherwise identical 
parameter values. 
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Figure 4. Directed transport in terms of individual scattering processes in parameter 
space, weak (panels a,c) vs. strong driving (b,d) for two-parameter, eqs. (|ll2p . (panels 
a,b) and single-parameter driving, eqs. (|1I3I4I) (c,d). Parameter values are Vq = 0.75, 
xo = 0.05 (a,c), 0.4 (b,d), uj = 1.0, ?A = 7r/2 (a,b), 7 = 0.1, <j) = tt/A (c,d). Color code: 
transmission 1 — >■ r and reflection r — > r blue, transmission r — )• 1 and reflection 1 — > 1 
red, else white. 



We analyze the resulting directed transport "locally in phase space" in figs. HI 
ini contrasting low with high driving amplitude (fig. H]), weak with strong symmetry 
breaking (l5H,b). and slow with fast driving (l5fc,d). Asymmetric scattering that 
contributes to directed currents, i.e., exit to the left (right), irrespective of the incoming 
direction, is labeled red (blue), neutral scattering (same outcome from both sides) white. 
The emerging structures refiect the relative amount of asymmetric scattering processes 
in terms of their total well as the complexity of the underlying phase space by 

forming fractals whiich resemble self-similar attraction basins [35] . 

Global evidence for directed transport is provided in fig. [6l comparing the 
dependence of the current on the incoming energy for various values of strength (panel 
a) and frequency u (b) of the driving. As a general tendency, it increases with both 
parameters but saturates as u approaches the inverse time of ballistic fiight through the 
scattering region, confirming our expectation that the mixed character of the dynamics 
is crucial for current generation. This is also consistent with the fact that currents 
vanish towards high energy: For Ei^ ^ Vq, trajectories belong to the regular part of 
phase space corresponding to quasi-free motion, cf. fig. El hence are always transmitted. 
At lower energies, current sets in as the first trajectories are refiected for only one of the 
incoming directions and reduces again as the same begins to occur also for the other 
direction. Superposed on these main trends, we observe frequent sign changes. The 
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Figure 5. Directed transport in terms of individual scattering processes in parameter 
space, weak vs. strong symmetry breaking (panels a,b) and slow vs. fast driving (c,d) 
for single-parameter driving, eqs. (jll3l4p . Panel a: weakly broken symmetry, 7 — 0.01, 
b; strongly broken symmetry, 7 — 1.0, c: low frequency, w = 0.05, d: high frequency, 
uj 10. Other parameter values are xq = 0.1, uj = 1.0 (a,b), 7 — 0.1 (c,d), </> — tt/2. 
Color code as in fig. |4l 
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Figure 6. Average currents as functions of various parameters. Single-parameter 
driving, eqs. (|l|3l4p . Panel a: current vs. initial energy for different driving amplitudes, 
xo — 0.01 (full red line), 0.1 (dashed, black), 0.5 (dotted, blue), 1.0 (dash-double 
dotted, green), 2.0 (dash-dotted, brown). Other parameter values are 7 = 0.1 x xq, 
(j) — 7r/2, UJ = 1.0. Panel b; current vs. initial energy for different driving frequencies 
Lo — 0.5 (full line, red), 1.0 (dashed, black), 5.0 (dotted, blue), 10.0 (dash-double- 
dotted, green). Other parameters are Vq — 0.75, a;o — 0.1, 7 — 0.1, (j) — '^1'^- 



corresponding sensitive dependence on the parameters [29] [32] allows for a fine tuning 
of the current. 

As is evident from figs. S] [5^,b, and [6^, even relatively weak driving or symmetry 
breaking suffices to achieve considerable transport. Conversely, this indicates that 
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currents due to small unintended asymmetries will be hard to avoid in a real laboratory 
setup. 



3. Quantizing chaotic pumps 

3.1. Floquet scattering theory 

Scattering at time-dependent potentials is generally inelastic so that standard quantum 
scattering theory no longer applies. If the driving is periodic, however, large parts of that 
familiar setting remain intact. A rigorous framework for the treatment of periodically 
time-dependent quantum systems is provided by Floquet theory [36l |37l [3H] , worked out 
for scattering systems in Refs. [391 SO] and applied to quantum pumps in [5]. We here 
only summarize a few facts directly relevant for us: 

• In periodically driven scattering, energy remains conserved mod/ia;, the photon 
energy associated to the periodic driving. For incoming and outgoing asymptotic 
energies this means 

-£^in/out(tt, "-in/out) = + "-in/out^, (6) 

where is the quasienergy pertaining to a Floquet eigenstate a and nin(out) 
enumerate incoming (outgoing) discrete Floquet channels. As a consequence, the 
total energy difference Ei = Eont — -Ein = l^ is quantized, I = Uout — '^in counting 
the total number of photons exchanged with the field. 

• The time evolution can be composed as a concatenation of discrete steps of duration 
T generated by the unitary Floquet operator 



Up = Texp 



dtH{t) 







(7) 



where H{t) is the Hamiltonian and T effects time ordering. They give rise to a 
discrete dynamical group and reduce numerical simulations of wavepacket scattering 
to repeated applications of Up [39l HO] . 

The "on- quasienergy- shell" scattering matrix •S'nj^^ ,^^^^(e) [391 HQ] is also organized 
in terms of Floquet channels Ui^, Uout- For one-dimensional scattering systems, it 
consists of four blocks S"'^ , a,T = — (+) denoting the left (right) asymptotic region. 
We define partial transmission and reflection probabilities as 

Tr'^iE.^+E,) = \S:^Z.JEin)\', nriE^n+E,) = \SZ,n.jEin)\M8) 
and the associated total quantities as 

oo oo 

T^riE^n) = Yl Tr^iEin), RZ{E,.) = Yl RriEin). (9) 

l=—oo l=—oo 

In terms of eqs. ([9]), the pumped probability current from left to right for fixed E 
is obtained in analogy to the classical probability flow ([5|) as 

HE,,,) = T-JiE;,,) + i?++(i?i„) - T+-(Ei,) - i?ro7(^in). (10) 
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If we assume the left and right leads to connect to reservoirs, the charge current is 
given as an average weighted by respective energy distributions f"{E) |5], 

{I)e = f E / [T~+{E)f-{E)-Tr{E)r{E+E,)].{ll) 

For Fermi distributed electrons at zero temperature and identical chemical 
potentials fi = Ep, this means 

{I)e = t E / [^r^(^) - T^iE + E,)]. (12) 



3.2. Numerical methods 

The matrix elements 5'^i^,nout (^) basic input to further evaluation of currents 
fll0|lllll2p etc. are determined numerically by wavepacket scattering: An initially 
free wavepacket l^/^J^) at a well-defined wavenumber /co(e) = \/2me/h is propagated 
stroboscopically through the scattering region by repeated application of f/p. Once a 
suitable termination criterion is met (e.g. when the occupation probability within the 
interaction region, cf. the integral on the r.h.s. of eq. (lT5l) . has decayed to below some 
threshold value), the final wavepackets \ipo\it) decomposed into Floquet channels at 
ki{Ein) = A/2m(i?in + Ei)/h, cf. eq. ([6]), to determine the scattering matrix elements 
^0,1 = ko{Ei^)i)l^^{ki{Ei^))/ki{E;^)ip?^{ko{Ei^)). The Floquet operator is efficiently 
calculated by the (t, t')-method [H]. 

Spurious interferences between the reflected and the transmitted wavepacket can 
occur if periodic boundary conditions are imposed at the ends of the finite spatial 
box underlying the propagation procedure. To prevent them, we employ absorbing 
boundaries instead. Specifically, the so-called Smooth Exterior Scaling Complex 
Absorbing Potentials (SES-CAPs [12]), while cumbersome to handle numerically, keep 
complications due to the corresponding non-Hermitian terms in the Hamiltonian at a 
minimum, as they are energy independent and leave the physical potential unsealed. 



3.3. Quantum transport 

In order to make contact with the adiabatic approach to pumping and to demonstrate the 
strongly nonadiabatic character of chaotic pumps, we compare in fig. [Th- numerical results 
for the two-parameter model ([2]) in the Floquet approach (full lines) with the predictions 
of the adiabatic scattering approach based on the Berry phase [3j (broken) : For a system 
with a single degree of freedom, two contacts (left, right), and two independently driven 
parameters ^i(t), ^2(^)5 the probability current through the device is 

hom = ^ I d'^BitE). (13) 



The area A in parameter space enclosed by the path ^(t) = (Ci(^), C2(^)) (cf. inset in 
fig. [6^) is penetrated by the "magnetic field" 

B{^, E)=\^2SZ{E)^^SZ\E) - ^^SZmE)^2SZ\E)] , (14) 
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Figure 7. Quantum current as a function of tlie initial energy. Panel a: Floquet 
approach, eq. (|T0| (full lines), compared with adiabatic theory, eqs. (|13ll4p (dashed), 
for slow (black, w = 0.1) vs. fast (red, lo = 1.0) two-parameter driving, eqs. (jll2p . 
Inset: parameter cycle in (a;_, a;+)-space (bold red line) and fictitious magnetic field 
B{x-,xj^), eq. ((T4)) at iJ = 0.45, underlying the adiabatic approximation (color code: 
from black, _B < 0, through red, B — Q, through yellow, B > Q). Panel b: Total current 
for single-parameter driving, eqs. (|ll3l4p . with uj — 1.0, — tt/2 (full line, red), Floquet 
currents in the elastic and first and second inelastic channels, n = 0, 1, 2 (dotted black, 
dot-dashed green, and dashed blue lines, resp., the latter two amplified by factors 10 
and 5 x 10^, resp.), and dwell time, eq. ([T5l) (short-dashed brown, rightmost ordinate). 
Other parameters are Vq = 0.75, xq = 0.1, h — 0.5. 




Figure 8. Fermi-averaged quantum current (IT2l) for single-parameter driving, 
eqs. (Ill3l4p . Panel a: Mean current vs. Fermi energy Ey for frequencies w = 0.1 
(fuU line, red), 0.3 (dashed green), 0.6 (dashed blue), 1.0 (dotted pink). Panel b: 
Scaled mean current I/uj vs. frequency uj at Fermi levels E-p — 0.2 (full line, red), 0.4 
(dashed green), 0.6 (dashed blue), 0.8 (dotted pink), 1.0 (dash-dotted bright blue). 
Other parameters areVb = 0.75, xq = 0.1, h = 0.5. 



with ^Sl^) = En,n' Sn]n'{E), and di = d/d^i. While for slow driving (black, w = 0.1), 
we find appreciable agreement, the results deviate drastically for fast driving (red, 
u = 1.0), indicating that strongly nonlinear transport mechanisms are involved. 

Adiabatic theory fails for single-parameter pumps, featured in fig. [Tb- We 
decompose the total current (full line) in Floquet channels (dashed). Although the 
main contribution comes from the elastic channel {n = 0), inelastic scattering [n = 1,2) 
is also present. Sharp resonances in the partial and total currents coincide largely with 
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peaks in the dwell time (dotted brown), measured as the fraction of the wavepacket 
inside the interaction region summed over stroboscopic time |28] . 



Sign changes are observed as well. The current even persists upon Fermi averaging (1121) 
(fig. [8]). As in the classical case (fig. |5]d), it increases with frequency as with energy but 
saturates for a; ~ 1 and E-p ^ \. Scaling the current as {I)e/<^ (fig- Eb) reveals this 
systematic deviation from the proportionality I ^ u, expected from linear response, for 

4. Conclusion 

Pumping in the non-adiabatic regime of fast and strong driving where nonlinear classical 
dynamics becomes relevant, remains largely unexplored. As a first survey into this realm, 
on the classical as well as on the quantum level, we have studied directed transport 
induced by irregular scattering. In particular, we have pointed out and provided 
evidence for the possibility of generating currents by driving just a single parameter, 
which manifestly goes beyond the adiabatic approximation. Their sensitive parameter 
dependence reflects the underlying chaotic dynamics and suggests to be harnished for 
control purposes. The classical nonlinearity induces similar nonadiabatic quantum 
transport: According to quantum adiabatic transport theory [3j, single-parameter 
current generation is impossible. Hence its occurrence for sufficiently strong and/or 
fast driving cannot be explained by the mere breaking of TRI alone. Analyzing this 
quantum-classical relation in the numerically very demanding regime of small effective 
Planck's constant where semiclassical methods apply |13] remains as a future task. 

Our model has been kept utterly simple, focused on the analysis of chaotic scattering 
and its impact on directed transport. Several options are conceivable to include 
more realistic details: Many-particle phenomena like charge blocking, dissipation, 
decoherence, and finite-temperature energy distributions in the reservoirs are obviously 
relevant. To be adequately treated, they require sophisticated nonequihbrium methods 
such as the Keldysh Green function technique [TOl HH US] . Other tempting perspectives 
to pursue are incorporating internal degrees of freedom, such as in particular spin |32] . 
as well as pumping against a potential gradient [26] and the rectification of noisy input 
forces. 
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